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Abstract 

The compact Abelian Higgs model is simulated on a cubic lattice where it possesses 
vortex lines and pointlike magnetic monopoles as topological defects. The focus of 
this high-precision Monte Carlo study is on the vortex network, which is investigated 
by means of percolation observables. In the region of the phase diagram where the 
Higgs and confinement phases are separated by a first-order transition, it is shown 
that the vortices percolate right at the phase boundary, and that the first-order 
nature of the transition is reflected by the network. In the crossover region, where 
the phase boundary ceases to be first order, the vortices are shown to still percolate. 
In contrast to other observables, the percolation observables show finite-size scaling. 
The exponents characterizing the critical behavior of the vortices in this region are 
shown to fall in the random percolation universality class. 

Key words: Compact Abelian Higgs Model, Monte Carlo Simulations, Vortex 
Network, Percolation 



1 Introduction 



The Abelian Higgs model with a compact gauge field formulated on a three- 
dimensional (3D) lattice possesses an intriguing phase structure [T|2|3] . In addition 
to the Higgs state where the photon acquires a mass, it exhibits a state in which 
electric charges are confined. The richness of the model, which serves as a toy 
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model for quark confinement, stems from the presence of two types of topological 
excitations, viz. vortex lines and magnetic monopoles. The latter are point defects 
in three dimensions which arise because of the compactness of the U(l) gauge 
group. In the pure 3D compact Abelian gauge theory, monopoles are known to 
form a plasma that physically causes confinement of electric charges for all values 
of the inverse gauge coupling (5 [1]. Being the only parameter present, the pure 
model therefore possesses only a confinement phase. The coupling to the scalar 
theory preserves this confinement state and gives in addition rise to a Higgs state. 
For sufficiently small values of the Higgs self-coupling parameter A, the two ground 
states are separated by a first-order transition [5|6] as sketched in Fig. [H For A 
larger than a critical value Ac, which depends on the value of the gauge coupling, 
the two states are no longer separated by a phase boundary across which thermo- 
dynamic observables become singular, as was first shown by Fradkin and Shenker 
[1] in the limit A — > oo where fiuctuations in the amplitude of the Higgs field be- 
come completely suppressed. In other words, it is always possible to cross over from 
one ground state to the other without encountering a thermodynamic singularity. 
Because of this, the Higgs and confinement states were thought to constitute a 
single phase, despite profound differences in physical properties. 

This conclusion is supported by symmetry arguments [7] . The relevant global sym- 
metry group of the compact 3D Abelian Higgs model (cAHM) is the cyclic group 
Zg of q elements, where the integer q denotes the electric charge of the Higgs field. 
For the doubly charged case (g = 2), the relevant symmetry group is Z2, which is 
in agreement with the known result that the model undergoes a continuous phase 
transition belonging to the 3D Ising universality class [T1I8] . For g = 1, this argu- 
ment excludes a phase transition characterized by a local order parameter in the 
spirit of Landau because the group Zi, which consists of only the unit element, 
cannot be spontaneously broken. Stated differently, there exist no local order pa- 
rameter which distinguishes the Higgs from the confinement state. 

In Ref. [H], we argued that the phase diagram is more refined than implied by 
this picture. We conjectured that although analytically connected, the two ground 
states can be considered as two distinct phases. The nature of the phase bound- 
ary is intimately connected to the distinct physical properties of the Higgs and 
confinement phase. For the latter phase to confine electric charges, the monopoles 
must form a plasma. This in turn can only happen when the line tension of the 
vortex lines, or fiux tubes, connecting monopoles and antimonopoles vanishes, so 
that they are no longer tightly bound in pairs as in the Higgs phase (for typi- 
cal configuration plots see the snapshots in Fig. [1]). Since the vortex line tension 
is finite in the Higgs phase and zero in the confinement phase, we argued that 
the phase boundary is uniquely defined by the vanishing of the vortex line ten- 
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Fig. 1. Sketch of the phase diagram of the Abelian Lattice Higgs Model as a function 
of the Higgs self-coupHng parameter A and the hopping parameter k at fixed inverse 
gauge couphng /3. For A < Ac the transition from the confined to the Higgs phase is 
of first-order. For A > Ac ordinary observables show no singular behavior and signal a 
crossover between the two ground states indicated by the dark region emanating from 
the point where the first-order transition line terminates. 

sion, irrespective of the order of the phase transition. The confinement mechanism 
operating in the 3D cAHM is essentially the dual superconductor scenario PUIITT] . 

In addition to open vortex lines, each having a monopole and an antimonopole at 
its endpoints, the system also possesses closed vortex lines. These are expected to 
be characterized by the same line tension as the open lines. Because of the finite 
line tension, large vortex loops are exponentially suppressed in the Higgs phase. 
Upon approaching the phase boundary, the line tension becomes smaller so that the 
vortex network can grow larger and the overall line density increases. Finally, at the 
phase boundary where the line tension vanishes, vortices can grow arbitrarily large 
at no energy cost. The phase boundary between the Higgs and confinement phase 
is therefore expected to be marked by a proliferation of (open and closed) vortex 
lines, as was first proposed by Einhorn and Savit [12]. The vortices proliferate both 
in the region where the transition is first-order and in the region where it is not. 
A line along which geometrical objects proliferate, yet thermodynamic quantities 
and other local gauge-invariant observables remain nonsingular has become known 
as a Kertesz line [13] . Such a line was first discussed in the context of spin clusters 
in the 2D Ising model in an applied magnetic field. 

The same deconfinement transition driven by the proliferation of vortices with a 
phase boundary consisting of a first-order line which ends in a critical point and 
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then continues as a Kertesz line was originally proposed by Langfeld for the SU(2) 
counterpart of the Abelian Higgs model defined on a 4D lattice [H] . The relevant 
vortices, forming surfaces in 4D, are center vortices which carry a flux related to the 
nontrivial element of the Z2 center of the SU(2) gauge group. The presence of the 
Kertesz line in this non-Abelian model has been further numerically investigated 
and confirmed in Ref. [T3]. While the concept of a Kertesz line was originally 
introduced in the context of lattice gauge theories in Ref. [TB], the interpretation 
of a deconfinement transition as driven by percolating vortices was previously put 
forward in the context of the SU(2) lattice Higgs model pT|18fl9] . and the pure 
SU(2) lattice gauge theory [201121122] . 

The purpose of this Monte Carlo study is to investigate the vortex proliferation 
scenario suggested in Refs. pi23] by studying the behavior of the vortex network 
directly. Because vortices are geometrical objects, their analysis is amenable to 
the methods developed in percolation theory [21]. We conjectured in Refs. [U|23] 
that along the Kertesz line, percolation observables have the usual percolation 
exponents. In addition, we expect that the vortex network displays discontinuous 
behavior in the region where the phase boundary consists of a first-order transition. 
For other recent Monte Carlo studies focusing on various aspects of the model see 
Refs. [^261127] and references therein. 

A similar high-precision Monte Carlo study of the behavior of a vortex network 
was recently carried out of the 3D XY and the lattice model [28f29] . respec- 
tively. The latter model, whose critical temperature on a cubic lattice is known to 
high precision, corresponds to taking /3 00 in the Higgs model so that the gauge 
fields become completely ordered. An important observation made in that study 
was that the overall vortex line density behaves similar to the energy, and the asso- 
ciated susceptibility similar to the specific heat. The vortex percolation threshold 
estimated through finite-size scaling analysis of the overall line density data was 
found to be perfectly consistent with the critical temperature. However, estimates 
based on any of the vortex percolation observables used, while being close to that 
temperature, never coincided with it within error bars. Possibly, the mismatches 
are related to the way the vortex networks are traced out, or to the absence of 
a stochastic element as present in, for example, the Fortuin-Kasteleyn definition 
of spin clusters in the Potts model [30]. In any case, since we use the same per- 
colation observables and the same (imperfect) rules to trace out vortex networks, 
we expect to find a similar mismatch, at least in the region where the transition 
is no longer discontinuous. To remind the reader of these qualifications, we will 
sometimes refer to the estimated percolation threshold as "apparent percolation 
threshold" . 
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The rest of the paper is organized as follows. In the next section, the observables 
used to investigate the model are introduced. In Sec. [SI the simulation methods as 
well as the numerical tools used to analyze the data are discussed. In that section, 
also a new tool to effectively collapse data gathered on lattices of different size 
is presented. Sections H] and [5] contain the Monte Carlo results in the vicinity of 
the phase boundary for the two regions where it consists of a first-order transition 
and a Kertesz line, respectively. In Sec. El the dependence of the location of the 
percolation threshold is investigated as a function of the parameters of the model. 
Finally, Sec. [7] contains a discussion of the Monte Carlo results and our conclusions. 



2 Definitions and Observables 



The Abelian lattice Higgs model with compact gauge field at the absolute zero of 
temperature is defined by the action S = Sg + S^f,, with the gauge part 

Ss = P [l-cosV(^)], (1) 

where P is the inverse gauge coupling parameter, P = with a the lattice spac- 

ing and q the electric charge of the Higgs field. We exclusively consider the case 
g = 1. The doubly charged Higgs field (g = 2), which has an even richer topolog- 
ical structure than the q = 1 case, has recently been investigated in Ref. [3T] . 
where it was found that the monopoles form chains. The sum in Eq. ([1]) ex- 
tends over all lattice sites x and lattice directions /i, and 6^y{x) denotes the 
usual plaquette variable 9^^{x) = A^9,^{x) — A^9^{x) with the lattice derivative 
A.u9^{x) = 9^{x + u) — 9fj,{x) and the compact link variable — tt < 9^{x) < tt. The 
matter part S*,^ of the lattice action is given by 

(2) 

where polar coordinates are chosen to represent the complex Higgs field = 
p{x)e^'^^^\ with the compact phase — vr < <f{x) < n. In Eq. (I2D, k is the hopping 
parameter, and A the Higgs self-coupling. We study the theory on a cubic lattice, 
which either is taken to represent a 3D space or spacetime box, depending on 
whether one of the dimensions of the lattice is interpreted as (Euclidean) time. 

In addition to measuring field observables such as the total action or energy S, the 
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hopping energy 

R,. = 

L3 



= 75 H p{x)p{x + fi) COS [A^(^(x) - , (3) 

X,fM 



the square Higgs amphtude, and coshnk energy 

^ = -T7^ E cos [A^^(x) - qe.ix)] , (4) 

we especially probe for topological excitations. A gauge invariant vortex line seg- 
ment jx{x) pointing in the A = 1, 2, 3 direction is given by 

jxix) = ex^vV^ly{x) - nx{x), (5) 

where l^{x) is the integer-valued field related to the electric current along the links 
of the lattice via 



= ^ {V^^(x) - qe^{x) - [V^<^(a;) - q6^,{x)]2n} , (6) 



and nx{x) measures the multiples of 271 with which the plaquette variable is shifted 
away from the interval [vr, vr): 

nxix) = -^(^x^iu {0^u{x) - [9^y{x)]2T,} ■ (7) 

Here, we use the usual modulo operation [a]27r = a — 2Txn which subtracts n 
multiples of 2tt from the variable a such that [a] takes values in the interval 
[7r,7r). While l^{x) measures the quantized vorticity, nx{x) gives the number of 
elementary Dirac strings piercing the plaquette with its normal pointing in the A 
direction. 

Monopoles are detected by taking the divergence of j^, m(x) = V^j^(x), where 
m{x) takes on integer values only. Using these definitions, we record the vortex 
line and monopole densities 

^ = ^EbM(^)U M^i^^m(x). (8) 

^ x,^i X 

As already mentioned in the Introduction, in the |0|^ theory, v behaves similar 
to the energy and the associated susceptibility similar to the specific heat [29]. A 
short summary of results for the observables not involving vortices can be found 
in Ref. 0. 
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Fig. 2. Networks formed by vortex lines. The small spheres represent monopoles and 
antimonopoles. Left: Higgs phase, characterized by a few short vortex lines which can 
either be closed or open. Right: Confinement phase, characterized by a large percolating 
vortex network. 

The main focus in this paper is on the vortex networks formed by individual vortex 
lines (see Fig. Tracing out such a network is ambiguous. We restrict ourselves 
to the simplest convention by defining a vortex line, which can be either open or 
closed, as a set of connected vortex line segments. Line segments are said to be 
connected if they enter or leave the same lattice cube. With this convention, four 
or six line segments entering and leaving a single cube are not further resolved 
into separate vortices, but are lumped together into one, self-intersecting vortex. 
The size tt, of a vortex is just the number of links forming the vortex. Vortex line 
segments with > 2 are not split into distinct segments and are only counted 
once in the vortex density. The vortex networks found with the help of these 
tracing rules are therefore only an approximation to the true networks. And the 
observed, or apparent percolation threshold will in general differ from the true 
one. We trace out the vortex networks by using a recursive algorithm and employ 
standard percolation observables to probe them [21]. Specifically, we determine 
the percolation probability Pid by recording when a vortex spans the lattice in an 
arbitrary direction. This definition is such that a vortex is already said to percolate 
when it spans the lattice in just one, arbitrary direction. We also determine the 
probability P^d that a vortex spans the lattice in all directions, and the percolation 
strength Poo, defined as the size of the largest vortex per lattice site. We have in 
addition considered the average vortex size, but this data is not well suited for 
analyses as we also found in previous studies on simpler models [32]. Finally, we 
use the susceptibility 
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(without the L for Pid and Psd) and the Binder parameter 



B{0) = 1 - ^ (10) 



of an observable O to probe the phase boundary. 



3 Simulation and Data Analysis Techniques 



3.1 Monte Carlo 



A variety of Monte Carlo methods are applied to efficiently simulate the system 
in the different parts of the phase diagram. In the first-order transition region, 
primarily the multicanonical algorithm (MUCA) [33] is used with a weight iteration 
as described in Ref. |34]. Weights are rendered flat in the hopping term ([3]) of the 
action to enhance tunnelling. In addition to local updates, the Higgs amplitude 
and the gauge angles are also updated globally [35] to allow for larger jumps in 
phase space and thus for shorter tunnelling times between the two metastable 
states. In the vicinity of the Kertesz line, the gauge fields are updated using the 
Metropolis algorithm, while the Higgs field is updated by means of heat-bath and 
overrelexation algorithms [36] . 

Measurements are typically taken after each sweep of the lattice and the entire time 
series is recorded to allow for error analysis and post-simulation data processing. 
Table [T] provides an overview of parameters used in the simulations. An integral 
part of our analysis are Ferrenberg-Swendsen reweighting techniques [37]I3S] . These 
techniques are applied in two distinct flavours. In the Kertesz region, we use the 
multihistogramm reweighting form, which means that we combine simulations at 



Table 1 

Number of sweeps of the lattice of linear size L and the number of different values of k 
considered for each lattice used in the simulations. 



/3 


A 


L 


sweeps 


different k values 


1.1 


0.2 
0.75 
0.025 


10, 12, 14, 16, 18, 20, 22, 24, 26, 30, 32, 42 
8,10,12,16,20 
8,10,12,14, 16,18,20 


10^ - 10^ 
10^ 
10^ 


20 
20 

1 (MUCA) 


2.0 


0.2 


8,10, 12,14,16,20,24,28,36 


10^ 


10 



8 



0.25 



0.15 




0.05 - 



0.473 0.4735 0.474 0.4745 0.475 0.4755 0.476 0.4765 



Fig. 3. Data points for the susceptibilities x(Ad) and x{C) (multiplied by 400) for 
P = 2.0 , A = 0.2 and different values of k and their interpolating lines obtained through 
reweighting. The measurements are taken on a cubic lattice of linear size L = 36. 



different parameters in an optimal way. Standard reweighting with MUCA weights 
is applied in the first-order transition region. To apply these techniques in the best 
possible way we use optimization routines such as the Brent method [39M0] to 
search for maxima in susceptibilities or crossing points of Binder parameters. The 
use of the Jackknife method [H] on top of these methods allows for error estimation 
of the results so determined. 

As an example of this approach, we show in Fig. [3] the susceptibilities of the 
percolation probability Pid and the coslink energy C. Note that the peak height 
of x(Pid) is exactly 1/4, since (Pf^) = (Pid) and hence x(-Pid) = (-Pid)(1 - 
(Pid)) is maximal at (Pid) = 1/2- The points in the plot correspond to individual 
simulations and the lines through the data points are obtained by reweighting. 
Different lines in the figure correspond to having different Jackknife blocks (but 
the same block in each time series) omitted. Notice that the peak in the percolation 
susceptibility is well defined whereas the peak in the coslink susceptibility is rather 
broad. The flatness of the peak is reflected in larger error bars on the estimated 
peak location. 

We have carefully checked our methods against known results for limiting cases 
such as the London limit A — >• oo, where fluctuations in the Higgs amplitude are 
completely frozen, and the XY model, obtained by setting in addition j3 = oo. 
A good check whether vortex lines and monopoles are correctly identified is to 
perform a gauge transformation under which their locations are invariant. 
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3.2 Finite-size scaling methods 



We estimate critical exponents through finite-size scahng analyses. The scaling 
Ansatz for an observable O states that in the vicinity of a continuous phase tran- 
sition, the data obtained for different values of the tuning parameter k and for 
different lattice sizes L fall onto a (weakly) universal scaling function Jq defined 
through 

OM=L'^/''fo{tL'/''), t={K-K,)/K,. (11) 

Here, denotes the critical point, v the correlation length exponent, and \o is 
the critical exponent characterizing the observable O. The variable t is the reduced 
coupling and measures the distance from the critical point. 

Data collapse is usually a good check whether the right critical exponents and 
critical couplings are found by other means. With the correct values, the measured 
data should fall on the universal curve given by Eq. (fTTI) . Here, we reverse this 
idea. Starting from an initial guess for the critical exponents and couplings, we 
compute the rescaled observable Ol{x) = Ol{k)L~^°/'' with x = tL^^" and judge 
the quality of the collapse in the interval [xmin, a^max] by introducing the weight 
function 

al= dx 01{x)-0l (x) , (12) 



with Ol{x) = J2l C^L{x)/nL and rii the number of different lattice sizes included 
in the sum perfect collapse in the window [xmin, a;max] would correspond to 

(T^ = 0, whereas a bad collapse has a large a^. This approach qualifies for imple- 
mentation as an automated method when combined with optimization algorithms 
to minimize by adjusting the values of the critical exponents. The method 
requires that data in each point in the interval [xmin, a;max] of the rescaled x-axis 
be compared, most of which has not been measured. A possibility is to interpo- 
late between data points by using a polynomial expansion of the function fa and 
to fit its coefficients to the data Ol [ISj. The fit then allows for calculating ap. 
We use a different approach in that we apply the standard reweighting techniques 
mentioned above to calculate Ol at every point k and hence Ol at every point 
X. We feel that this approach is more natural and does not add further degrees of 
freedom to the process. 

In practice our implementation is as follows. We first reweight the measured data on 
all system sizes. We then start with initial values for the exponents and couplings, 
and apply a minimization algorithm, such as the simplex method [39113], that 



varies the exponents until is minimized. Since the algorithm can become locked 
in a local minimum, the procedure must be repeated for many starting points. We 
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find that for our purposes, where only three parameters need to be determined, the 
routine gives rehable and consistent results. To obtain error estimates, we apply 
the Jackknife method on top of the whole process. The effect of correction terms to 
scaling is minimized by repeating the procedure for increasingly smaller intervals 
[xmin, a^max] arouud the point where collapse is attempted. We expect our tool to be 
useful in other simulation studies as well. An idea similar to our implementation 
was recently presented in Ref. |44j . 



4 First-Order Transition Region 

We start our numerical study by investigating the behavior of the vortex network 
in the first-order transition region of the phase diagram sketched in Fig. [H We 
choose to simulate at /5 = 1.1 and A = 0.025 where the first-order transition 
is strong enough already on small lattices so that time consuming simulations 
on larger systems are not needed. The left plot in Fig. H] shows time series of 
the hopping energy and the percolation strength normalized by the volume 
as measured in a typical MUCA simulation [33] . Changes in the energy, which 
reflect the more or less random walk through the metastable region, are seen 
to occur jointly with changes in the percolation strength. The right plot shows 
the correlation histogram of the number of configurations without a percolating 
network versus hopping energy E^. Larger negative energies are seen to strongly 
correlate with the absence of a percolating vortex network while smaller negative 




75000 76000 77000 78000 79000 80000 -2 -1.8 -1.6 -1.4 -1.2 -1 -0.8 -0.6 



MC time Ei, 

Fig. 4. Left: Time series of a MUCA simulation showing tiie hopping energy E^^ (upper 
curve) and percolation strength Pqo (lower data points) for L = 16. The two observables 
are seen to be intimately connected. To fit in one plot, the hopping energy is shifted 
upwards. Right: Correlation histogram of the number of configurations without a per- 
colating network over hopping energy E'h for system sizes L = 10 and 16. Note the 
logarithmic scale used on the vertical axis. 
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Fig. 5. Left: Binder parameters B{Ey^) (upper curves) and B{Poo) (lower curves) as a 
function of k in the first-order transition region at (3 = 1.1, A = 0.025 measured on 
lattices of linear size L = 10,12,14,16,20 (from broad to sharp peaks). The Binder 
parameter B{Poq) has been divided by 100 so that both parameters fit into one plot. 
Right: Scaling plots of the maxima of the various Binder parameters and susceptibilities. 
The Binder parameters show larger finite-size corrections. All curves converge within 



error bars to the same infinite-volume value of k 



per 



0.424570(3). 



energies strongly correlate with the presence of such a network. 

The central thesis we put forward in Ref. [9] is that the phase diagram of the model 
can be understood in terms of proliferating vortices. In the first-order region, the 
location of the phase transition has been determined to high precision with the 
help of observables not involving vortices such as the hopping energy. To determine 
the location of the percolation threshold in the infinite- volume limit, we consider 
the vortex percolation probability Pid and percolation strength Poo, and analyze 
the scaling of the locations Kper{L) of the maxima in the associated susceptibilities 



and Binder parameters with lattice size L. At a first-order transition, 
expected to scale as 



"per 



"per 



Oil 



(L) is 



(13) 



with c a constant. We reweight the MUCA time series to obtain the susceptibilities 
and Binder parameters in the vicinity of K-per{L), and use the methods presented 
above to search for the peak locations and its error bars. Figure [5] summarizes our 
results. It shows the scaling of the percolation thresholds K-per{L) obtained from 
the percolation probability Pid and the percolation strength Pqo together with the 
critical points Hc{L) obtained from the hopping energy i?h and the monopole den- 
sity M. The K-per{Ly^ obtained from the network observables are seen to strongly 
depend on the system size, while the fi;c(L)'s obtained from x{^) ^^e almost in- 
dependent of volume. The results for x(-^h) are not included in Fig. [5] as they 
cannot be distinguished from those for M on the scale of the figure. For lattice 
sizes larger than L = 12, all curves converge to the same point within error bars. 



12 



This brings us to the important conclusion that the vortex percolation thresh- 
old is located precisely at the first-order transition point. Our best estimate for 
the transition point, based on both percolation and nonpercolation observables, 
is Kc = 0.424570(3) which is determined from the average of the fit results sum- 
marized in Table O Averaging the estimates for percolation and non-percolation 
observables separately, we obtain Kper = 0.424572(4) and Kc = 0.424565(4), respec- 
tively. The absence of a mismatch between the critical temperature and the vortex 
percolation threshold, as was found in the |0|^ theory using the same percolation 
observables [29], is probably because the phase transition is discontinuous here. An 
abrupt change in the ground state apparently washes out any inaccuracy caused 
by an imperfect tracing out of the vortex network. 

To establish in an unbiased fashion that the vortex network reflects the first- 
order nature of the transition aX f3 = 1.1 and A = 0.025, we assume nothing 
about the transition and apply standard finite-size scaling to the vortex percolation 
probability Pid as for a continuous transition. According to Eq. fll3p . we then 
expect to find u^cr = 1/d = 1/3. To obtain an estimate for the percolation threshold 
Kper on the infinite lattice and Upcr, we reweight the data points from the MUG A 
time series and subsequently apply our collapsing routine with Kpor and z/por as 
free parameters. Figure E] shows the input and the result of this procedure. The 
extracted values, Upcr = 0.424568(6), t'per = 0.325(10), are perfectly consistent 
with a first-order transition precisely at the expected location. For comparison. 
Fig. [6] also shows the hopping energy which is seen to display the same behavior 

Table 2 

Results of fitting the data in Fig. [5] to Eq. ()13p . The fits starting at L = 8 and L = 10 
also include the correction term Those starting at L = 12 include only the leading 

term except for B(Poo), where the correction term is still needed to obtain reason- 

able results. The number in square brackets denotes P^r degree of freedom (DOF), 
X^/DOF. Results are grouped according to percolation and non-percolation observables. 



Observable 


L>8 


L > 10 


L > 12 


x{Pid) 


0.424569(5) [1.6] 


0.424569(5) [2.0] 


0.424570(6) 


[2.3] 


XiPoo) 


0.424572(6) [1.9] 


0.424572(5) [2.2] 


0.424573(6) 


[2.5] 


B{Poo) 


0.424608(10) [10.7] 


0.424608(10) [13.1] 


0.424577(6) 


[1.1] 


X{M) 


0.424567(3) [2.4] 


0.424566(5) [3.1] 


0.424567(6) 


[4.2] 


x{Eh) 


0.424564(5) [2.52] 


0.424564(6) [3.1] 


0.424564(6) 


[3.2] 


B{M) 


0.424558(5) [2.6] 


0.424558(8) [3.5] 


0.424555(5) 


[3.2] 




0.424565(6) [3.3] 


0.424568(5) [3.9] 


0.424566(6) 


[4.0] 
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Fig. 6. Left: Reweighted hopping energy £'h and percolation probability Pid for various 
lattice sizes as a function of the hopping parameter close to the first-order phase tran- 
sition. The crossing points should provide a second method to estimate the transition 
point. The estimates obtained using the two observables agree within error bars. Right: 
Collapse of the Pid data with Kpcr = 0.424568, f, 
Eh data (not shown) is achieved for Kc = 0.424563 and ly = 0.338. 



per — u.-i^iuuu, lyper — 0.325. An Optimal collapse of the 



as the percolation probability. Repeating this procedure for the hopping energy 
data, we obtain = 0.424563(6) and u = 0.338(10). Both these estimates based 
on the data collapse analysis are perfectly consistent with the previous estimates 
from finite-size scaling. The raw data in Fig. [H] suggest as in [13] that the crossing 
points of the curves measured on lattices of different size provide a second method 
to estimate the transition point. 



5 Kertesz Line 



We continue our analysis of the vortex network in the region where the transi- 
tion ceases to be of first order. In Ref. [9], we postulated that in this part of the 
phase diagram the Higgs and confinement phases are separated by a Kertesz line. 
Along this line vortices proliferate, yet thermodynamic quantities remain nonsin- 
gular across it. The conjecture is based on the numerical observation that in the 
Higgs phase, the monopoles are tightly bound in monopole-antimonopole pairs. 
The magnetic flux emanating from a monopole is squeezed into a magnetic flux 
tube (vortex) which ends on an antimonopole. The flnite line tension forces the 
vortex lines to be short. In the conflnement phase, the monopoles are no longer 
bound in pairs, but form a plasma. For this to arise, the vortex line tension must 
vanish. Vortex lines, both open and closed, can then grow arbitrarily long at no 
energy cost and proliferate. To facilitate comparison with our previous work [9], 
we choose the parameters (3 = 1.1 and A = 0.2. Figure [7] (left) shows our results 
for the vortex line density as a function of the hopping parameter k for lattices 
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of linear size varying from L = 12 to L = 26. The inset gives the corresponding 
susceptibihties, which display the same remarkable behavior first observed in the 
London limit A — oo for other observables [8p5] . Namely, the susceptibility data 
obtained on lattices of different sizes is seen to collapse without rescaling. In par- 
ticular, the maxima of the susceptibilities do not show any finite-size scaling. From 
these data we obtain the estimate = 0.5615(7) for the location of the phase 
boundary. Applying the same analysis to the hopping energy, we arrive at the esti- 
mate Kx = 0.5655(7) which is close, but does not agree within error bars with the 
previous estimate. To further investigate this issue, we consider the monopole den- 
sity to see if it leads to the same estimate as does the hopping energy. The results 
of our initial study indicated that both estimates agree within error bars. How- 
ever, the high-precision Monte Carlo simulations carried out in the present work 
show that this is not the case for the monopole density yields Kx = 0.5639(1). In 
evaluating these results, it should be kept in mind that none of these susceptibil- 
ities diverge in the infinite- volume limit, so that small discrepancies were to be 
expected. The change in the ground state is for this reason usually referred to as 
a crossover (x). The discrepancy between the peak locations of these observables 
is shown below in Fig. O 



In Ref. [9], we argued that the phase diagram is more refined than just show- 
ing a crossover between the Higgs and confinement ground states in that a sharp 
boundary between the two phases does exist in the form of a Kertesz line across 
which the vortices proliferate. Moreover, as we argued partly on the basis of sym- 
metry, the percolation observables in the vicinity of the Kertesz line should be 
characterized by the usual percolation exponents. Unlike the observables previ- 
ously studied, these observables are expected to show finite-size scaling. Figure [7] 
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Fig. 7. Left: Vortex line density and susceptibility (inset) as a function of the hopping 
parameter k for systems of linear size L = 12,16,18,26, and /? = 1.1, A = 0.2. Right: 
The susceptibility x(-foo) as a function of the hopping parameter k for systems of linear 
size L = 22,26,32,42. 
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(right) shows the susceptibihty x(-foo) of the percolation strength as a function of 
the hopping parameter k for lattices of linear size varying from L = 22 to L = 42. 
It is indeed observed that this percolation observable depends on the lattice size 
even though we considered system sizes larger than those for which the other 
observables already reached the infinite- volume limit. That is, in contrast to the 
other observables, percolation observables allow for a precise location of the phase 
boundary. 

To estimate the percolation exponents, we study the behavior of the percolation 
probability Pid and P3D as well as the percolation strength Pqo in the vicinity of 
the percolation threshold /tper (see Fig. [H]). On the infinite lattice, the percolation 
strength vanishes on approaching the threshold as Poo ~ (k — Kper)^'"'% while 
x(Poo) diverges as x(-Poo) ~ |k — ^pcr | . Finally, the correlation length ^per, 
which provides a typical length scale of the vortex network, diverges as ^per ~ 
— Kpcrl"'^'""'- Given the discrepancy found in the context of the |0|^ theory, we do 
not expect the estimate of the percolation threshold using percolation observables 
to coincide with the one based on the vortex line density. We therefore determine 
the location of the percolation threshold anew together with the exponent v^er by 
studying the finite-size behavior of the susceptibilities of percolation observables. 
Figure [9] shows the scaling of the locations of the susceptibility maxima with 1/L. 
The observables considered are the percolation probabilities Pid and P3D, and the 
percolation strength P^o- For comparison, also the data for the hopping energy, 
which was used in Ref. [9] to estimate the location of the phase boundary, and the 
data for the monopole and vortex densities are included. The percolation data in 
Fig. [9] are fitted to the function 

fi:per(L) = Kper + cL-i/'^p-, (14) 



using the standard least-squares method. 




Unfortunately, this approach does not 
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Fig. 8. Percolation probability Pid and percolation strength as a function of the 
hopping parameter k for systems of linear size L = 22, 26, 30, 42 at /3 = 1.1, A = 0.2. 
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Fig. 9. Locations of the susceptibility maxima of the percolation probabilities Pid and 
P3D, and the percolation strength Pqo as a function of l/L for /3 = 1.1. The dashed line 
is a fit of Xmax(-foo) to Eq. (I14p . The horizontal line is through the corresponding data 
for the vortex density, while the two sets of data points above this line pertain to the 
monopole density and hopping energy. 

give reliable and consistent values for Kper and fper when repeated for the differ- 
ent observables and lattice sizes. The most stable results are obtained from the 
susceptibility of Poo, giving Kper = 0.54141(18) and z/pcr = 0.87(6). Results for Pid 
and P3D depend too much on the fitting regime and we conclude that Eq. ffT^ for 
these observables is only fulfilled for large lattice sizes. As expected, the estimate 
of the location of the percolation threshold does not agree with Kx obtained from 
the vortex line density. 

To improve these estimates, we now apply our collapsing routine, rather than 
carrying out additional simulations for other lattice sizes. The method has the 
advantage that much more data is used as input for not only the locations of the 
peak maxima but all the data in the vicinity as well as the interpolated values 
obtained from reweighting are included. We take the Binder parameter B{Poo) of 
the percolation strength together with the scaling Ansatz 

P(Poo) = h {tL^/'^--) , t={K- Kpe,)/Kper, (15) 

where his a, scaling function. Figure [10] shows the raw data together with the data 
collapse for k^^v = 0.54145 and z/per = 0.881. The good quality of the data collapse 
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0.54145(2) and i/per = 0.881(2). 



is apparent from the 3D plot in the same figure. Given these estimates for Kper and 
'per, we next apply the collapsing routine to determine the exponents /5per and 7per 



from the scaling relations 

Poo = 



X 



(Poo) = (tL^/"-^ 



(16) 



with / and g scaling functions. Table [3] summarizes our results and compares them 
with the random percolation exponents. In agreement with our conjecture [S], the 
estimates, which as shown in Fig. [TT] lead to a good data collapse, are consistent 
with standard percolation exponents. Notice that, in contrast to the findings in 
the |(/)|^ theory, the vortices appear to proliferate in the confinement phase after 
passing through the crossover region {k^ct < i^x)- We expect this to be related to 



Table 3 

Critical exponents of the percolating vortex network across the Kertesz line for f3 
and A = 0.2 compared with the standard percolation exponents. 
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Model 


'^pcr 




/3pcr 


Tpor 
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0.54145(2) 


0.881(2) 


0.43(2) 


1.76(2) 


Percolation [l6] (see also |47|48j) 




0.8765(16) 


0.4522(8) 


1.7933(85) 
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Fig. 11. Left: Collapse of the percolation strength Poo yielding an exponent 
/3p^^ = 0.43(2). Right: The same for the susceptibility x,{^oo) yielding ^per — 

1.76(2). 

the presence of monopoles which can act as sources for the vortices. To assess this, 
we study in the next section how the relative positions of the percolation threshold 
and the crossover region vary with changes in the parameters A and f3. 



6 Fine Tuning 

Since the vertex percolation threshold Kper and the phase boundary Hc coincide in 
the first-order transition region, it is expected that as we approach this region, the 
discrepancy found in the previous section at f3 = 1.1 and A = 0.2 becomes smaller. 
To verify this, we repeat our analysis at /3 = 1.1 and A = 0.075, which is closer to 
the first-order transition region. We find that the difference between the locations 
of the percolation threshold and the crossover determined using the overall vortex 
line susceptibility x(f) indeed becomes smaller, changing from about 0.018 at 
A = 0.2 to 0.007 here. Moving in the opposite direction of increasing A, we observe 
the discrepancy to increase, becoming as large as about 0.1 in the London limit 
A ^ oo. The percolation threshold still appears in the confinement phase after 
passing through the crossover region. That is, changes in the Higgs self-coupling 
A seem to leave the relative positions of the crossover and the apparent vortex 
proliferation threshold unchanged. 

We next vary P and set P = 2.0 and A = 0.2. By increasing P, one suppresses 
the monopoles. They completely disappear in the limit P oo, where the theory 
looses its compactness. Unfortunately, simulations at P = 2.0, A = 0.2 are com- 
putationally much more challenging than at /5 = 1.1,A = 0.2 as autocorrelation 
times are much longer. We therefore restrict ourselves to lattice sizes up to L = 36. 
Moreover, observables other than percolation observables become less useful as can 
be seen from our example in Fig. El Whereas the susceptibility of the percolation 



19 



0.4760 
0.4755 
0.4750 
0.4745 
0.4740 
0.4735 
0.4730 
0.4725 
0.4720 
0.4715 



- X X X X ^ 

---i; * * 


s 

X 


-43 


"it-. I 




- 


x{Eh) 






X{C) 






X{v) - ^ 













0.02 



0.04 



0.06 

1/L 



0.08 



0.10 



Fig. 12. Locations of the susceptibility maxima of various observables as a function of 
1/L for p = 2.0. 

probability has a pronounced peak and a clear maximum, the maximum in the 
coslink susceptibility would be difficult to identify without reweighting. The 
large error bars obtained for the peak location of x{C) reflect the absence of a 
pronounced peak. 



Our main conclusion of the simulations at /3 = 2.0 and A = 0.2 is that the loca- 
tions of the apparent percolation threshold and the crossover determined using the 
overall vortex line susceptibility x(^) have changed relative positions. This con- 
clusion is based on Fig. [12] in which Kper{L), estimated using different percolation 
observables, is plotted as a function of 1/L to see their tendency for L — > oo. 
The first observation is that the data points are much closer to each other than 
was the case at (3 = 1.1 (see Fig. As before, Hpcr{L) obtained from the sus- 
ceptibility of the percolation strength (lowest set of data points) shows the largest 
corrections but increases monotonically for L ^ oo. The values Kper{L) obtained 
from the percolation probability Pid show less drastic corrections but are more 
difficult to extrapolate to the infinite- volume limit. The two sets of data points 
both extrapolate to a value around Kp^r ~ 0.4748 which follows from ffiting x{Poo) 
to Eq. (HI). 

Figure [12] displays in addition to the locations of the susceptibility maxima of 
the vortex line density, also those of the coslink energy and the hopping energy. 
For L > 16 the heights and locations of these susceptibilities remain constant, so 
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that, as far as these observables are concerned, the infinite- volume hmit is reached. 
The location of the crossover in the infinite-volume limit determined using these 
observables is below all estimates of the percolation threshold. In other words, 
while for (3 = 1.1, A = 0.2, Upcr < Kx, here the relative positions have changed, 
Kper > Kx, and the vortices proliferate already in the Higgs phase before enter- 
ing the crossover region. By continuity we then expect at some value of (3 in the 
interval 1.1 < (3 < 2.0 the percolation threshold of the vortices to coincide with 
the location of the crossover. Since varying j3 physically changes the monopole 
density, it is tempting to conclude that monopoles impede the formation of per- 
colating vortex lines, and that by adjusting the monopole density, the location of 
the apparent vortex percolation threshold can be fine-tuned to coincide with that 
of the crossover. The use of the word "apparent" here is to underscore that vortex 
networks identified with our tracing rules are only an approximation to the true 
networks. The element of impediment introduced by the monopoles possibly plays 
a similar role as the stochastic element in the Fortuin-Kasteleyn construction of 
spin clusters in the Potts model |30j . 

We finally estimate the critical exponent z/per for j3 = 2.0. Using the percolation 
probability as input observable for our collapsing routine, we arrive at the value 
^pcr = 0.88(1) and Kpcr = 0.4748(1). As a crosscheck we fit the scaling of x(-Poo) 
with the lattice size L to Eq. 1^ giving Up^, = 0.88(8) and Kp^r = 0.47484(15) 
with a x^/DOF = 2.1. This is again consistent with the standard percolation 
exponent. Although the relative positions of the apparent percolation threshold 
and the crossover have changed, variations in the inverse gauge coupling parameter 
f3 appear not to change the value of this critical exponent. 



7 Conclusions 

The vortices arising in the compact Abelian Higgs model have been investigated 
by means of Monte Carlo simulations on a cubic lattice, and analyzed with the 
help of observables known from percolation theory. Because their behavior is more 
pronounced, percolation observables are better suited than other observables to 
probe the phase boundary, both in the first-order transition region as well as in 
the crossover region of the phase diagram. In the region where the Higgs and con- 
finement phases are separated by a first-order transition, the vortices percolate 
right at the phase boundary. Since the rules applied to trace out the vortices re- 
sult in networks that are in general only an approximation to the true ones, it is 
concluded that the discontinous first-order transition is forgivable of the resulting 
inaccuracies. The vortex network refiects the first-order nature of the transition in 
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this region of the phase diagram. In the crossover region, the vortices still percolate. 
The percolation observables show second-order critical behavior along the Kertesz 
line that is characterized by the usual percolation exponents. The location of the 
vortex percolation threshold estimated using percolation observables does not co- 
incide with that of the crossover estimated using the vortex line density. Ideally, 
one would expect both to coincide. Also the |0|^ theory in 3D, which undergoes a 
continuous phase transition, shows a similar behavior [22] • Whereas the estimate 
based on the line density coincides with the critical temperature, the estimates 
based on any of the percolation observables considered do not agree within error 
bars. This discrepancy may arise because the vortex networks are not correctly 
traced out or because a stochastic or impeding element in the construction of a 
network is missing. The monopoles appear to play such a role. 
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